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The nature of the critical point of the Anderson transi- 
tion in high magnetic fields is discussed with an emphasis 
on scale invariance and universality of the critical exponent. 
Special attention is paid to the distribution function of the 
conductance which becomes size and model independent at 
the critical point. The fractal properties of the wave function 
which are related to scale invariance are also discussed. 



I. INTRODUCTION 

The Anderson transition (AT) has been attracting the 
attention of condensed matter physicists for more than 
four decades (!],(§. The AT is a zero temperature quan- 
tum phase transition separating metallic and insulating 
phases which is induced by a spatially fluctuating ran- 
dom potential. The transition can be described using 
the scaling theory of localization [|||4|] . Near the critical 
point in three dimensional (3D) systems behavior which 
is typical of quantum phase transitions is observed for 
quantities such as the conductance g, correlation length 
£ etc. For example, as we approach the critical point by 
changing a parameter w such as strength of disorder or 
Fermi energy, the correlation length diverges as 



W - W c 



(1) 



while the conductivity a vanishes from the metallic side 
according to the power law, 



| to - w c \ s 



(2) 



If we approach the transition from the insulating side, 
then the dielectric constant e diverges as 

e~ \w~w c \- s '. (3) 

As in the critical phenomena of magnetic systems, the 
exponents u, s and s' are not independent but are related 



= (d - 2)v, s' = 2v. 



(4) 



A knowledge of v is enough to fix the critical exponents 
s and s' . These latter exponents can be measured exper- 
imentally though there has been controversy concerning 
the correct values pj . It is thus important to have a pre- 
cise theoretical estimated of v in order to compare with 
the experiments. 



The value of v is expected to be universal, i.e., inde- 
pendent of the details of the model and dependent only 
on basic symmetries such as that under the operation of 
time reversal. The classification of the critical behavior 
according to the symmetries of the system was predicted 
from field theoretic considerations [QU , and recently ver- 
ified numerically fl. 

At the Anderson transition the correlation length di- 
verges and the wave function becomes scale invariant. 
This invariance is characterized by a fractal dimension 
D%. This is reflected in the size independence of the dis- 
tribution function of the conductance HllOl as well as the 
statistics of the energy levels 
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In this paper we report a numerical simulation of the 
wave function dynamics at the critical point for the tight 
binding model in a magnetic field. We have observed 
anomalous diffusion and estimated the fractal dimension- 
ality. The universality of the distribution function of the 
conductance at the critical point has also been verified. 
The final section is devoted to summary of our results 
and concluding remarks in connection with experiments. 



II. ANDERSON TRANSITION IN MAGNETIC 
FIELDS 

Magnetic fields have two effects on the Anderson tran- 
sition. One is to delocalize the electronic states by break- 
ing time reversal symmetry, and the other is to localize 
the electronic states by shrinking the wave function due 
to cyclotron motion. Which of these two dominates de- 
pends on the situation. 

The tight binding Hamiltonian which incorporates the 
effect of magnetic fields is given by 



H = V 



E 

<i,j> 



(5) 



where C\{Ci) denotes a creation (annihilation) operator 
of an electron at the site i. Energies {Wi\ are distributed 
independently and uniformly in the range [-W/2, W/2]. 
By fixing the Fermi energy E to be, e.g., E = while 
increasing W, the system is driven to be an insulator at 
W = W c where W c is the critical disorder. The Peierls 
phase factors exp(i#ij) describe magnetic fields. The 
hopping amplitude V is assumed to be the energy unit, 
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V = 1. We assume a simple cubic 3D lattice for simplic- 
ity, and all the length scales are measured in units of its 
lattice constant a. 

In the absence of magnetic fields, the AT occurs at 
a critical disorder W c ~ 16.5 at the center of the band 
E = 0. The critical exponent v has been estimated as 
v = 1.59±0.03 0. Applying strong magnetic fields to the 
tight binding model it has been shown that the value of 
the critical point as well as the scaling curve change. In a 
magnetic field the value of v is estimated to be 1.43 ±0.04 

H . The exponent is not dependent of the strength of the 
magnetic fields || and is unchanged in random magnetic 
fields which can be realized by assuming a random phase 
for the hopping elemnts [fijif . 

Another model which describes the 3D Anderson tran- 
sition in high magnetic field is a stack of two dimen- 
sional layers where the strong quantizing field is applied 
perpendicular to the plane [r^| . In purely two dimen- 
sional systems, the quantum Hall effect (QHE) occurs. 
The electronic states are delocalized only at the centre 
of each Landau band. In this case the critical exponent 
^qhe is estimated to be close to 7/3 |L8]]. The introduc- 
tion of interlayer hopping between the layers makes the 
delocalized region finite in energy, and changes the ex- 
ponent to 1.45±0.15 |l!|. This is close to the value of 

I. 43 given above in agreement with universality of the 
AT. Universality has also been verified for changes in the 
transfer integral between the layers [ fl7j ]. 

There is another interesting aspect of the stacked layer 
model. If we impose Dirichlet boundary condition in- 
stead of periodic boundary conditions, we find mag- 
netic edge states circulating along the perimeter of the 
two dimensional layer. Stacking the layers along the z- 
direction, we then have novel electronic states where the 
electron can hop both in the +z and — z directions but 
can rotate only in the clock-wise or anticlock-wise direc- 
tion (which is determined by the direction of the field) 
in the x-y plane. The edge states which compose this 
strange "sheath" are critical That is to say, as 

we increase the linear dimension of the plane, say L, the 
localization length in the z-direction £|| also increases in 
proportion to L and diverges in the thermodynamic limit. 



III. CRITICAL BEHAVIOR 

A. Anomalous diffusion and fractal dimensionality 

In the metallic regime, the electron diffuses and the 
mean squared diffusion radius r 2 {t) = (t Ir 2 ]t) is propor- 
tional to the time t, 



- 2 {t) = 2dDt, 



(6) 



where D is the diffusion coefficient and the average is 
taken over disorder. In the insulating regime where the 



wave function ip are localized as exp(— r/£), the squared 
diffusion radius saturates at [Elf 



d(d + 1) ^ 
km r 2 {t) = £ , 

t—>oo 4 



(7) 



To understand the intermediate region we use the renor- 
malization group, from which 

7HT) =b 2 f((w-w c )b 1 /",tb- z ) 

where b is the scale factor in the renormalization group 
and z is the dynamical exponent. Time t is measured in 
units of h/V. From this equation we deduce the scaling 
form El 



r 2 (t) = t 2/z F(t 1/zu (w-w c )) 



(8) 



A similar relation holds for classical percolation theory 
p3| . On the metallic side of tke transition we expect that 
at sufficiently long times a linear in t growth of the mean 
squared radius in accordance with (||). For this to be so 
we must have F(x) ~ x s when x 3> 1 with s — (z — 2)v. 
Thus on the metallic side of the transition at long times 
we have 



r 2 (t) ~ \w - w c \ s t 

Since according to the Einstein relation a ~ D we see 
that s is indeed the exponent in (Q). For non- interacting 
electrons z — d and we recover the Wegner scaling law 
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{d-2)v 



On the insulating side of the transition we expect (0) to 
hold at long times. Imposing this in (ph leads to 



r 2 (t) - \w-w c \- 2 " 

confirming that v in (||) is indeed the exponent governing 
the divergence of the localization length. Exactly at the 
critical point w — w c and we see that the square diffusion 
length grows as 



(9) 



Another quantity which can be used to investigate the 
dynamics of wave function at the critical point is the 
return probability 
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C(t) = - / dt'\{t' 



t 



(10) 



which is related to the fractal dimensionality of the wave 
function D2 as 



c(t) ~ r D *i d . 



(ii) 



Direct diagonalization of the 3D systems requires huge 
CPU power, especially when the Hamiltonian is complex 
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because of the applied magnetic fields. Instead, we have 
used the equation of motion method to study the diffu- 
sion process. We prepare a wave packet |0) with the en- 
ergy E located at the center of the system, and calculate 
numerically the time evolution using \t) = e~ lHt ^ n \0). 
When evaluating the factor e - lHt / h ^ we use the decom- 
position formula for exponential operators |2^j2^]. In 
figure 1, we show the results for r 2 (t) calculated for a 
59 x 59 x 59 cubic lattice. The magnetic field is parallel 
to the z-direction and the magnitude of the flux per unit 
cell is 0.1 times the flux quantum. The critical disorder in 
this case is W c = 17.8 We see clearly the t 2 / 3 -law for 
r 2 (t) , confirming that validity of the scaling equation (||) 
as well as the scaling relation equation ([|) . The estimate 
of L>2 from C(t) is shown in the inset of Figure 1. We 
find Z?2 = 1-7 which is significantly smaller than the spa- 
tial dimension 3, demonstrating that the wave function 
at the transition is not at all similar to a typical extended 
wave function. This value is consistent with the recent 
estimate of Z?2 for layered systems in high magnetic field 



B. Conductance distribution 

In a d-dimensional hyper cubic lattice the dimension- 
less conductance g is defined as 



G ah 



d-2 



9 = 



(12) 



e 2 /h e 2 /h 

From the Landauer formula |28|] , we have 

g = trtt f (13) 

where t is the transmission matrix. The matrix t is ob- 
tained by iteratively calculating the Green function. [^9| . 

In 3D metallic LxLxL systems the conductance distri- 
bution function P(g) is the normal distribution, the mean 
of which (g) is proportional to the size L. The variance, 
on the other hand, is universal, a phenomenon which is 
known as universal conductance fluctuations f30|| . In the 
insulating regime P(g) is log-normal. 

At the Anderson transition, not only the variance 
but the distribution function itself becomes universal 
|p|, ^0pl]l . The resulting distribution function is plotted in 
figure 2 for L = 8, 10, 12 and 14. We plot P(log.g) instead 
of P{g) to see more clearly the detail of the distribution 
function. The histogram is for the uniform magnetic field 
while the dots are for the random phase hopping model. 
The critical disorder W c depends on the strength of the 
field and how we break the time reversal symmetry (i.e., 
uniform magnetic field or random phase hopping). Nev- 
ertheless, P(logg) at the critical points is universal. 

Once the system is away from the critical point, the 
P(g) begins to show size dependence. To demonstrate 
this size dependence, we plot P(log<?) in vanishing field 



away from criticality at W = 17.5 in figure 3. This value 
of the disorder is about 6% larger than W C {B = 0) = 
16.5. We can see subtle but clear size dependence of 

P(g)- 

Similar behavior of the conductance distribution is also 
observed in the layered system in high magnetic fields, 
though in this model, the system is highly anisotropic 
and the form of P c (g) is different ]32]. 



IV. SUMMARY AND CONCLUDING REMARK 

In this paper we have discussed several features of 
the Anderson transition which are related to the self- 
similarity of the eigenstates at the critical point. The 
fractal dimensionality of these critical eigenstates is al- 
most half the original space dimension, « 1.7. The square 
diffusion length r 2 (t) has been shown to grow as i 2 / 3 in 
the non-interacting model, irrespectively of the values of 
the critical exponents. At the transition, the distribution 
function of the conductance P(g) becomes model and size 
independent. The distribution function of g close to AT 
had been recently obtained experimentally which is 
consistent with our results. Such universality of the dis- 
tribution is also seen in the statistics of the energy levels 
at the transition [p7T|~p~5|| . 

We have also seen that the layered system in perpen- 
dicular magnetic fields shows interesting transport prop- 
erties. The critical behavior of this system may still be 
the same as that found in the tight binding isotropic sys- 
tem. 

In order to relate our results with experiments at finite 
temperature T, we now discuss the T dependence of the 
conductivity, cr(T). At finite temperature, the inelastic 
scattering time T; n and the inelastic scattering length 2; n 
are finite. In the metallic regime, they are related by l{ n ~ 

1/2 

r in . At criticality, the diffusion is anomalous diffusion 
and this becomes 

k-r,* (14) 
The effective diffusion coefficient D e s becomes 
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P>eS T- 



(15) 



leading to the conductivity at finite temperature <j(T) 

a(T) ~ (16) 

where we have assumed r m ~ 1/T. Setting s = v gives 
a ~ T 1 / 3 which is independent of the values of the ex- 
ponents s and v. In the presence of the electron-electron 
interaction, the relation (^J) may no longer be valid and 
instead 



s= (d-2 -6)v 



(17) 
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should be used p4| . This leads to the suggestion that 
the exponent of the temperature dependence is different 
when the time reversal symmetry is broken. In experi- 
ments 1 35 the 1/3-power is widely observed. This means 
that even in the presence of electron-electron interaction, 
the relation s = v may not be modified significantly. 
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FIG. 2. Distribution function of the logarithm of the con- 
ductance g at the critical point. Triangles (A), diamonds (O), 
squares (□) and circles (o) corresponds to L = 8, 10, 12 and 
14 for the random phase hopping model, respectively. The 
histogram is for a system in uniform magnetic fields. 



FIG. 3. Distribution function of the logarithm of the con- 
ductance g when the strength of disorder is slightly larger than 
the critical value W c . Triangles (A), diamonds (O), squares 
(□) and circles (o) again corresponds to L = 8, 10, 12 and 14, 
respectively. 
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